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An accurate risk assessment for fatigue damage is of vital disk are discussed, 
importance for the design and service of today's turboma- 
chinery components. We present an approach for quantifying 

the probability of crack initiation due to surface driven low- INTRODUCTION 

cycle fatigue (LCF). This approach is based on the theory of Due to the necessity for a more flexible service of gas 
failure-time processes and takes inhomogeneous stress fields turbines, low-cycle fatigue (LCF) design has become of es- 
and size effects into account. The method has been imple- sential importance in today's gas turbine engineering. The 
mented as a finite-element postprocessor which uses quadra- end of LCF-life of an engineering component is often defined 
ture formulae of higher order. Results of applying this new by the initiation of a crack of a certain size. LCF is mostly 
approach to an example case of a gas-turbine compressor surface driven so that an LCF crack initiates on the compo- 
nent's surface. In the case of polycrystalline metal the grain 
structure has a great influence on the LCF failure mechanism 
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cycles to crack initiation, even under lab conditions. 

Standard design approaches often derive a predicted 
component life with respect to LCF from the average life 
times of the most loaded points on the component plus safety 
factors which account for the scatter band, size effects and 
uncertainties in the stress and temperature fields. However, 
this method can have a lack of sufficient accuracy which can 
result in designs that are too conservative or too optimistic, 
for example. In particular, the Coffin-Manson-Basquin equa- 
tion and Wohler curves play an essential role in reliability 
estimations regarding fatigue. Moreover, size effects have 
a significant influence on fatigue life and are mostly con- 
sidered by safety factors. For more detailed discussions of 
fatigue confer 0, 0, and 0. 

In this work we focus on LCF in conjunction with poly- 
crystalline metals and in particular on the number of load 
cycles until crack initiation. We present a local and proba- 
bilistic model for LCF according to and which we use 
to estimate the risk for LCF crack initiation on a compres- 
sor disk. Here, the reaction of a component to cyclic loads 
is taken into account via a linear elastic finite element anal- 
ysis (FEA) and via Neuber shakedown, see 0, and ||9l , 
respectively. 

Because LCF cracks are small in the initiation phase 
they will only influence stress fields on micro- and 
mesoscales and so it is reasonable to assume that crack for- 
mation in one region of the component's surface is not influ- 
enced by the crack forming process on another part. Thus 
crack formation can be considered as a problem of spatial 
statistics iflOll . Then, 13 and infer that hazard rates for 
crack initiation have to be integrals over some local func- 
tion depending on the local stress or strain fields. Therefore, 
hazard rates for the component can be expressed by a sur- 
face integral over some crack formation intensity function 
depending on local fields. The latter can also be regarded as 
the density for the intensity measure of a Poisson point pro- 
cess (PPP) in time and space, confer ifTTl and fPHl . In the 
role of the PPP is also emphasized in that context. 

Following and we model the hazard rates with a 
rather conservative approach by means of intensity measures 
of Weibull type. Note that scale-shape distributions are very 
common in reliability statistics, confer 0131 . Furthermore, 
we assume that the scale variable Adet is of the same func- 
tional form as the usual Coffin-Manson-Basquin equation, 
confer 0, and below. This results in a Weibull distri- 
bution for the number N of cycles of first crack initiation on 
the component as well. 

In contrast to this purely phenomenological approach 
avoids detailed modeling at the meso scale which facilitates 
calibration with experiments. Both approaches have the use 
of Weibull distributions in common. From a materials engi- 
neering point of view the model of and has the signif- 



icant advantage - compared to standard methods in fatigue 
- of bypassing the standard specimen approach and consid- 
ering size effects. Besides LCF-test results with standard- 
ized specimens, different strain-controlled LCF results can 
be used together for the calibration of the model such as re- 
sults from specimens with different geometries or under dif- 
ferent inhomogeneous strain fields. 

In order to numerically compute the Weibull distribu- 
tion, quadrature formulae are employed for the correspond- 
ing integration. As locations of stress concentrations result 
in higher nonlinearities in the integrand we use quadrature 
formulae of higher orders. For this purpose we interpolate 
the field values according to principles of mesh and finite- 
element generation, confer 0, lfl4l and lfT5l . 

Having obtained the Weibull shape and scale param- 
eter from calibration and from numerical integration, re- 
spectively, the corresponding distribution function yields the 
probability for LCF crack initiation with respect to the num- 
ber N of load cycles. In the design process one can decide 
which number N is acceptable corresponding to the assigned 
risk. Another important quantity of the LCF model is the 
local crack initiation density. This field shows how much 
each region of the surface contributes to the overall expected 
number of crack initiations and leads to critical as well as to 
possibly overengineered parts of the component. 

This paper combines some aspects of materials engi- 
neering, reliability statistics and FEA. In the first section we 
consider linear elasticity, Neuber shakedown, fatigue analy- 
sis and we present the local and probabilistic model for LCF. 
In Section [2] we first focus on FEA. Then, we discuss nu- 
merical integration and important functions of our presented 
approach. The last section shows results of applying our 
method to an example case of a gas-turbine compressor disk. 

1 A LOCAL AND PROBABILISTIC MODEL FOR 
LCF 

In this section we discuss results of linear isotropic elas- 
ticity, Neuber shakedown and fatigue analysis. The last sub- 
section presents the local and probabilistic model for LCF 
which is introduced and motivated in and 0. 

1.1 Linear Isotropic Elasticity and Neuber Shakedown 

In this work we consider designs made from single- 
phased polycrystalline metal. Usually metallic material does 
not consist of one single crystal but of different crystalline 
regions which are called grains. These grains have an or- 
der of magnitude typically in the micrometer and millimeter 
scale. The continuum mechanical approach assumes that the 
considered sizing scale is large compared to the inter-atom 
distances. Thus the material is considered to be smeared and 
all quantities are continuous. In general, single crystals have 



anisotropic properties. In polycrystals the orientation of the 
grains is randomly distributed. If the grains are sufficiently 
small compared to the component's size, the anisotropic ef- 
fects of the grains average out and an approximately isotropic 
material can be assumed. 

In the following, we employ the continuum mechanical 
approach and assume isotropic material behavior at scales 
significantly larger than the grain size and assume suffi- 
ciently small deformations. Thus, linear isotropic elasticity 
can be applied to describe the behavior of components from 
single-phased polycrystalline metal under external loading 
and plasticity can be considered by Neuber shakedown. In 
this section we present theoretical backgrounds and closely 
follow Section 2 of which is based on lfl6l and (T). Note 
that our example case of a compressor disk is subject to a 
homogeneous temperature field so that we do not consider 
thermoelasticity. 

Let £2 be a domain which represents the component 
shape filled with a deformable medium such as polycrys- 
talline metal which is initially at equilibrium. Moreover, let v 
be the normal on the surface dQ. of £2, let f be an external load 
and let u be the three-dimensional displacement field in £1. 
Finally, let dQ.D,dQ.N be a partition of the boundary where 
dQ.o is clamped and on dQ.N a normal load g is imposed. 
Then, according to [8 | the mixed boundary value problem 
(BVP) of linear isotropic elasticity is described by: 



V-o e (u) + f=0 in£2 



(1) 



with O e (u) = X(V ■ u)I + /j(Vu + Vu r ) and with boundary 
conditions u = on 3£2d and G e (u) • V = g on dQ.N- Here, 
X and jj are the Lame coefficients. The linearized strain 
rate tensor e e (u) is defined as e e (u) = j(Vu + Vu r ), i.e. 

efj = j (j^ + f° r hj = 1,2,3. Numerical solutions of 
the BVP can be computed by an FEA, confer fl6l . and 
Section|2]below. 

The knowledge of the threshold between elastic and 
plastic deformations is very important as plastic deforma- 
tions can allude to an imminent residual fracture. According 
to U this threshold is often described by so-called yield cri- 
teria. In this work we use the von Mises yield criterion which 
is given by 



[(<Ji - o 2 ) 2 + (oi - o 3 ) 2 + (a 2 - o 3 ) 2 ] = kf , (2) 



where kp is the critical value of the criterion and 01,02,03 
are the principal stresses. Here, the left-hand side is pro- 
portional to the elastic strain energy of distortion. If the 
criterion is applied to uniaxial tensile tests the relationship 



kf = R p / y/3 is obtained, where R p is the critical value of the 
only nonzero principle stress in a uniaxial tensile test. The 
von Mises stress is defined as 



[(oi - o 2 ) 2 + (01 - 03) 2 + (02 - 03) 2 ]. (3) 



Thus, the previous criterion can be written as o v = R p and be 
used to predict yielding of metal under any loading condition 
from results of uniaxial tensile tests. 

In the following we introduce the Ramberg-Osgood 
equation, confer ifPTl . It can be used to locally derive strain 
levels from scalar comparison stresses such as the von Mises 
stress. These strain levels determine strain-controlled fatigue 
life. The Ramberg-Osgood equation establishes stress-strain 
curves of metals near their yield points. It is very accu- 
rate in the case of smooth elastic -plastic transitions which 
can be observed for metals that harden with plastic defor- 
mations, for example. If K denotes the cyclic strain harden- 
ing coefficient and n the cyclic strain hardening exponent the 
Ramberg-Osgood equation is given by 



CTv 
E 



L/n 



(4) 



with Young' modulus E = >i \^ 2t '' . The equation defines the 
comparison strain £,,, where we also write £ v = RO(o v ). 

Finally, we present the method of Neuber shakedown, 
confer lfl"8l and JTJ . Let o e v denote the von Mises stress 
which is obtained only from linear elastic computations and 
let G y be the von Mises stress which also considers plastic- 
ity. Here, o v is called elastic-plastic von Mises stress. If lin- 
ear elasticity leads to stress values greater than the material 
yield strength, Neuber shakedown will estimate correspond- 
ing elastic-plastic stress values. An energy-conservation 
ansatz is the foundation of Neuber's approach. It results in a 
relationship between the elastic von Mises stres|3 <5 e v and the 
elastic -plastic von Mises stress o v : 



o v e v 



h O v — 

E \K 



l/n 



(5) 



Here, the Ramberg-Osgood approach is also used. K t is the 
notch factor, which is set to one if o e v is obtained from the 
BVP (03 where notches are incorporated in the boundary def- 
inition. Given the elastic comparison stress o^„ we can thus 
calculate the elastic -plastic von Mises stress by solving (0. 
Thereby, we are able to obtain £,, from (3). Note that we also 
write o v =SD- l (o e v ). 



1 Confer [19] for details on Neuber shakedown in conjunction with equiv- 
alent stresses. As an alternative to this method one could also use Glinka's 
method, confer e.g. (9). 
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Fig. 1 . EN-DIAGRAM OF A STANDARDIZED SPECIMEN. 

1.2 Fatigue and the Coffln-Manson-Basquin Equation 

Fatigue describes the damage or failure of material un- 
der cyclic loading, confer JT], ||20) and 0. Major examples 
of cases where fatigue occurs are activation and deactivation 
operations, e.g. of motor vehicles and of gas turbines, and 
oscillations in technical units. Material science analyzes the 
physical nature of fatigue and ways to determine the num- 
ber of cycles until a material fails under cyclic loading. In 
this work, we will consider a compressor component of a 
gas turbine subject to surface driven low-cycle fatigue (LCF). 
For backgrounds on surface driven LCF failure mechanism 
with respect to polycrystalline metal we refer to JT], ||2], 
and 0. 

We now consider important methods of fatigue analysis 
and closely follow Section 3 of Q and |UJ. In fatigue speci- 
men testing the number of cycles until failure is determined. 
If the tests are strain controlled so-called E —N diagrams - 
see Figure [JJ - are created, where the relationship between 
the strain amplitude e a and number iVj of cycles until crack 
initiation is called Wohler curve. Usually, the range of cycles 
is subdivided into low-cycle fatigue (LCF) and high-cycle fa- 
tigue (HCF). LCF loads are often strain controlled, whereas 
HCF loads are mainly stress controlled so that corresponding 
S — N diagrams are analyzed. 

For the purpose of analysis the strain amplitude z a is 
subdivided into an elastic and plastic part where e„ = £ e J + £„' 
holds. In the LCF range the plastic part £§' dominates, 
whereas in the HCF range the elastic part plays a greater 
role. Introducing the parameters fatigue strength and fa- 
tigue strength exponent b we present the so-called Basquin 

equation e e J = -^-(2A^,) / ' for the elastic part, where E is 
Young's modulus. Regarding the LCF range the Coffin- 
Manson equation e£ = e'f(2Ni) c describes the denominating 
plastic part, where the parameters and c are called fatigue 
ductility and fatigue ductility exponent, respectively. For a 
detailed discussion of the physical origin of this equation we 
refer to j4). The combination of the previous equations leads 



to the Coffin-Manson-Basquin (CMB) equation 

e a = e e J+E? 1 = -^{Wif + ^iWiY. (6) 

The parameters can be calibrated according to the test data - 
see FigureQ]- by means of maximum likelihood methods, for 
example. Confer l2T1 . |fl3l , |20l and Subsection |3.1| below. 

Structural design concepts with respect to LCF often 
consider the component's surface position of highest stress 
and then analyze the Wohler curve which corresponds to the 
conditions at that surface position. Mostly safety factors are 
imposed to account for the stochastic nature of fatigue, size 
effects^ and for uncertainties in the stress and temperature 
fields. Note that sometimes several surface positions of high- 
est stress are considered which depends on the component. 
This concept is called safe-life approach to fatigue design 
and is often used in engineering as well as very similar meth- 
ods, confer UJ. 

Note that several extensions exist to the CMB equation 
©. In particular, when approaching the HCF region, mean 
stress effects are of increasing importance. The modified 
Morrow equation is one approach that would consider such 
effects. For further discussions confer |0. 

1.3 From Reliability Statistics to Probabilistic LCF 

Now, we introduce the local and probabilistic model for 
LCF as presented in J5) and J6) which can be derived from 
reliability statistic, confer lTT3l . 

We model failure-time processes on continuous scale, 
although time in our context is a number of load cycles 
and thereby an integer number. Let N denote a contin- 
uous random variable which represents the time of crack 
initiation here identified with failure of a system or com- 
ponent. If P denotes the underlying probability measure 
^jv(w) = P(N < n) is the cumulative distribution function and 
/iv( n ) = dFf/(n)/dn the density function. The survival func- 
tion is defined by 5V(«) = P(N > n) = 1 — Fn(ji) and the 
hazard function by 

,, s ,. P(n<N <n + An\N>n) f N (n) 

h(n) = lim = 

v An^O An 1 - F N (n) 

confer (13\ . h is also called hazard rate or instantaneous fail- 
ure rate function. For a small step An the expression h(n) - An 
is an approximation for the propensity of an object or sys- 
tem to fail in the next time step An, given survival to time 
n. For a large number m(n) of items in operation at time n 



2 Note that different geometries of test specimens lead to different 
Wohler curves, confer (2]- 



the product m(n) ■ h(n) is approximately the number of fail- 
ures per unit time. Defining the cumulative hazard function 
H(n) = Jo h(t)dt one can show that the survival function sat- 
isfies Sff(n) = 1 —F}f(n) = exp(— H(n)). This shows that a 
model ansatz for the hazard function leads to a corresponding 
distribution function F^. 

Now, we introduce the crucial assumption for the local 
and probabilistic model for LCF of (5) in the case of poly- 
crystalline metal. Consider LCF failure mechanism on the 
component which is represented by the domain Q.. We as- 
sume that the surface zone that is affected from the crack ini- 
tiation process of a single LCF crack is small with respect to 
the surface of the component. This surface zone corresponds 
to faces of a few grains. As long range order phenomena are 
unusual in polycrystalline metal, we pass to the following 
assumption: 
Assumption (L) 

In any surface region A C 3£2 the corresponding hazard rate 
Iia is a local functional of the elastic displacement field u in 
that particular region with 

h A (n)= f p(n;Vu,V 2 uWA. (7) 

Ja 

Here, Vu is the Jacobian matrix and V 2 u the Hessian of 
u. Note that the loads that we consider in the given con- 
text (mainly elastic plastic stresses and strains) can all be ex- 
pressed as functions of Vu. The integrand p is called hazard 
density function. 

Assumption (L) is also called the property of spatial 
additivity of hazard rates. A motivation of (L) and more 
detailed backgrounds can be found in J5]. For mathemati- 
cal simplicity, we restrict ourselves to the case where only 
the dependence on elastic strains (or equivalently stresses) is 
taken into account. 

In the case of inhomogeneous strain fields assumption 
(L) implies h(n) = J sa p(n;E e )dA for some hazard density 
function p. Because of /*j\r(n) = 1 — exp(— H(n)) = 1 — 
exp(— /q h(t)dt) we obtain for the probability of failure in 
dQ. until cycle n: 

F N (n) = 1 -exp^- £ J^p(t-e e )dAd?j . (8) 

The ansatz ([8]) can also be derived by the Poisson point pro- 
cess with p(n;£ e (x)) as the intensity measure. For details on 
point processes we refer to ifTTI and lfl"2l . The advantage of 
the point process is that also the probability of a given num- 
ber of cracks initiations in A C 3£2 within n load cycles can 



be computed via 

_ z q 

P(number of crack initiations on A = q) = e — (9) 

q. 

for z = Jo J A p(t;s e )dAdt. But if cracks have grown suffi- 
ciently large they will mutually influence their local stress 
fields and thus the approach will break down. 

We now establish a link to deterministic LCF analysis 
via CMB equation (O which leads to an appropriate choice 
for the hazard density function p. We assume that the number 
N of cycles to crack initiation are Weibull distributed which 
can be realized by the choice of a Weibull hazard ansatz 

/ \ m— 1 

p(n;x) = p(n;e e (x)) = — —— ( ; — rrr 1 

(10) 

Here, m is the Weibull shape and A^j e t(e e ) the Weibull scale 
parameter which is supposed to depend on the elastic strain 
tensor e e (x)of the BVP |[T). Combining ^ and ( fTOt leads to 
the following model of : 

(Local and Probabilistic Model for LCF) 
Let the scale field A^j e t( x ) = -^det(£a( x )X x £ d£2, be the so- 
lution of the CMB equation © 

e a (x) - -^(2^V det (x)) i, + e}(2^V det (x)) c -, (11) 

where E a (x) is computed from e e (x) via^ linear isotropic 
elasticity, from the von Mises stress G v (x), from a fl (x) = 
SD~ (a v (x)/2) according to Neuber shakedown and from 
the Ramberg-Osgood equation with e„(x) = RO(o a (x)). 
Then, the local and probabilistic model for LCF is given by 
the cumulative distribution function 

Fjy(n) = l-exp(- p [ dAds). (12) 

for n > and some m > 1, which yields the probability for 
LCF crack initiation in the interval [0,n]. 

The shape parameter m determines the scatter of the dis- 
tribution where small values for m > 1 corresponcQto a large 
scatter and where the limit m — > °° is the deterministic limit. 
Note that the Weibull hazard function can be easily replaced 
by any other differentiable hazard function with scale param- 
eter A^ det . 



3 Confer Sections [TTlandfni 

4 < m < 1 is not realistic for fatigue. 



The CMB parameters of the model are not the same as 
obtained from fitting standard specimen data. We calibrate 
them by means of usual maximum likelihood methods, con- 
fer Section [5] ETl and lfl3l . Furthermore, note that volume 
driven fatigue could be considered as well by replacing the 
surface integral in (fTZt with a volume integral whose inte- 
grand only differs by different material parameters. For a 
discussion of volume driven fatigue such as HCF confer JT]. 

With respect to materials engineering, the local and 
probabilistic model for LCF has significant advantages com- 
pared to the safe-life approach to fatigue design: The model 
bypasses the standard specimen approach and takes size ef- 
fects into account, i.e. results from arbitrary geometries un- 
der LCF failure mechanism can be employed to calibrate 
our model and every position of the surface of an engineer- 
ing part is considered by a surface integral which does not 
need information on Wohler curves of a specific specimen. 
Thereby inhomogeneous stress fields are taken into account. 

2 FINITE ELEMENT ANALYSIS AND POSTPRO- 
CESSING 

In this section we describe our finite-element postpro- 
cessor which computes the distribution function with respect 
to fatigue life. First, we briefly introduce into concepts 
of FEA. Then, we discuss numerical integration and back- 
grounds on the postprocessor. 

2.1 Finite Element Analysis and Lagrange Elements 

In order to numerically solve the B VP ([T) of linear elas- 
ticity on a three-dimensional polyhedron £2 we apply FEA, 
where a so-called weak formulation is considered on a finite- 
dimensional space of functions, confer (8), lfl5l and lfl4l . 
This function space and the geometry £2 are described by a 
mesh of finite elements, where each element consists of a 
three-dimensional compact and connected set T and of a fi- 
nite set of basis functions n = {\\t\ \\t„ sh } associated to a 
set of nodes {ai, . . . ,a„ s/! } in T. In the following {7\n} de- 
notes a finite element and the functions of n are called shape 
functions. The FEA solution of the BVP ([T]i restricted to T 
is then a certain linear combination of the shape functions. 

Very popular finite elements are the Lagrange finite el- 
ements which are contained in most FEA packages. In our 
example case of a compressor disk we will use Abaqus 6.9-2 
and Lagrange elements of Serendipity class (C3D20) where 
the shape functions are of the form 

V|/(x)= £ «ii,/ 2 ,!'3-^i4 2 4 3 ( 13 ) 

0<i 1 ,i 2 ,/ 3 <2 
'l+i2+'3<3 



Now, principles of mesh and finite element generation 
are briefly introduced, where we closely follow . 
(Mesh) 

The mesh of a domain £2 is given by the union of com- 
pact and connected sets {T m }i< m <N el with £2 = Um=i T m 
and all interiors nonempty and pairwise disjoint, hj de- 
noting the greatest distance of two points in T we define 
"Kb = {T m }i<m<N d where h = m<at{h T | T G {T m }i< m <N d }- 

The starting point for the construction of meshes is a 
geometric reference cell T from which the different sets are 
generated. 

(Generated Mesh, Reference Cell) 
A mesh HQ, of a domain £2 is called generated if there exists a 
diffeomorphisirQ Tt for every T <G HQ, an d a fixed set f such 
that Tt(T) = T. The set T is called reference cell. 

Now, we define the so-called geometric transformation 
which can provide the previous diffeomorphism, confer : 

(Geometric Transformation) 
Let {r,n} be a finite element with shape functions 
and let {ai, ...,a„ , } be three-dimensional 
nodes. Then the geometric transformation denotes the map 

>hh 

S^T(S) = £a,%(«). (14) 

i=i 

In most cases a mesh generator produces a list of « geo 
times N e i nodes such as {af , — ,a^ geo }i< m <Ar ei . Here, N e i is 
the number of elements again. The vectors {a™, ...,3.™^} 
are also defined as the geometric nodes of the m-th element. 
If Y m is the geometric transformation corresponding to the 
reference cell T and to {a'", ...,&'"} with m = 1,... ,N e i we 
set T m = Y m (T). So we obtain m elements and can write 
down the following definition: 

(Geometric Reference Finite Element) 
If a finite element {T^FLec,} is used for mesh generation by 
means of given nodes 

{a'",...,a^ o }i< m < Nel , (15) 

its geometric transformations T„, and setting T m = Y m (T) for 
m= 1 , . . . ,N e / , it is called the geometric reference finite ele- 
ment. 

(Generation of Lagrange Elements) 

Finally, we consider the generation of Lagrange elements 
{r,FI} from a reference Lagrange element {r,n}. Here, 
{a,'}i<i<„. denotes the Lagrange nodes of {7\n} and 
. . . , \j/„ j;i } the corresponding shape functions which are 



for coefficients oc ;i ./ 2 ,i 3 determined by \|/y(ay) = 8y. 



5 A diffeomorphisms is continuously differentiable map whose inverse 
exists and is continuously differentiable as well. 



assigned to the nodes via %(ay) = 8,/. In this case fl = ft geo 
and given nodes {a™,...,a% }i<m<N e i provide a mesh % 

with reference cell T, geometric transformations Yt and 
sets T = Yt(T). The nodes of T satisfy a, = Ty(a,) for 
i = l,...,n s h- In order to construct II consider for functions 
v the map 



vi — ► <j)r (v) = v o Yt (16) 



which is linear and invertible. Then, the shape functions are 
given by n = {^ 1 {p)\p eft}. 



3D Finite Element Constructions: Tetrahedrons and Bricks 
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Fig. 2. GEOMETRIC TRANSFORMATION AND INTEGRATION 
POINTS ON THE UNIT TRIANGLE AND RECTANGLE. 



2.2 Postprocessing 

Recall that the local and probabilistic model for LCF is 
given by the cumulative distribution function (fT2] i. In the 
following we address numerical integration and explain how 
the linearized strain field £ c (x),x 6 £2, is computed by means 
of linear elastic results from FEA with Lagrange elements. 

The surface integral in (fTZt is numerically computed 
by quadrature formulae. Note that we can also numerically 
compute volume integrals, where in (fT~2T > d£l is replaced by 
£2. This can be of higher interest if we extent our model 
to volume driven failure mechanism such as HCF Because 
locations of stress concentrations will have a major contribu- 
tion to our probability functional we have to consider higher 
nonlinearities in the integrand by our numerical approach. 
Therefore, we will use quadrature formulae of higher order. 

Since £2 is a polyhedron which consists of finite ele- 
ments we can conduct the integration of (17\ on every face 
of a finite element which contributes to the boundary of dQ. 
and then sum up each of that integral values. Moreover, the 
geometric transformation (fl4l can be used as a chart to trans- 
form each face integration on the corresponding map area 
which is the unit triangle in case of tetrahedrons or the unit 
rectangle in case of bricks. Also see Figure [2] where exem- 
plary geometric transformations and integration points are 
depicted. For more details on surface integration consider 
the comments to (l22l below and the term of integration over 
submanifolds as described in 11221 . for example. Substitution 
for multiple variables is employed in case of volume inte- 
grals and thus the element integrals are transformed to the 
unit tetrahedron or unit brick. 

(Quadrature of Order k) 

Let K be an integration domain, let the l q numbers C0i , . . . , CO; 
be denoted as weights and the l q points £ K as 

integration points (quadrature points). The weights and in- 
tegration points are called quadrature of order k, if k is the 



largest integer such that 



f '" 
p(x)dx=T(i)ip(^i) for all 

Jk z t1 



a h >d x \ ■■■ x d 



°<h <rf<* 

h+-+id<k 



(17) 



«/; i d real 



With respect to intervals [a,b] Table [TJ shows corresponding 
weights and integration points. Confer Table 8.1. and 8.2. in 
(H) as well, where Table 8.2 contains quadrature formulae for 
triangles. Note that a quadrature on an interval can be used 
to obtain a quadrature on a rectangle [a,b] x [c,d] by sub- 
dividing the multidimensional integral into one-dimensional 
integrals. This results into the four integration points of the 
unit rectangle in Figure [2] Similarly we obtain quadrature 
formulae for bricks and thus can numerically integrate over 
volumes of finite elements. For volume quadratures on tetra- 
hedrons and quadratures of higher orders we again refer to 
Section 8 in GO. 

Considering the finite element generation, we use the 
coordinates of the nodes, the displacements of the nodes and 
the connectivity (information on which nodes belong to a 
specific element). Moreover, we have to employ the explicit 
form of the shape functions . . . , } of the reference 
Lagrange element {T,Tl}. Then, we obtain the geometric 
transformation for every finite element {7), II/} according to 
(fl4l i and the local displacements on every 7] according to 
(US: 



xi — >-(ur ( oT r; )(i) = £u (fe %(x), (18) 

k=] 

with ut displacement vectors at each node of {7), 11/}. 



Table 1 . QUADRATURE ON THE INTERVALS = 

(a + b)/2 AND 8 = b-a. 



\a,b] WITH to: 



1 1 



3 2 



5 3 



7 4 



Int. Points 



TO± 



2\'3 



m± 



m±§ v / (15+2V30)/35 



m±§ ( V(15- 2^30)735 



Weights CO; 




Recalling ( fT2l , we explain the computation of 
Let %> be the y'-th face of the i-th element for i = I,... ,N e i 
and j — 1,...,Nf, where Nf is the number of faces of the 
considered element type. Let S^-..^ be 1 if Jlj is a subset of 
the surface dQ. with area greater than 0, otherwise let it be 0. 
Then, 



mis 
JdnN dei \N det 

N e , N F r 

'=!;= i L 



m— 1 



rfAdfc = «'" / ) c/A 

Jan V^det , 



1 



dA 



(21) 



In case of bricks we integrate over the unit rectangle as map 
area: 



Now, the computation of the linearized strain rate tensor 
E e is explained. We will express the fields in dependence of 
x e T. Let x be the coordinates that describe the locations 
of 7}. Because of the geometric transformation we have x = 
T7j(x) for some xef. Recall the form of the linearized 
strain rate tensor: 



dxi 



Ue {1,2,3}, 



with u i components of displacement vector u. In order to ex- 
press the strain completely in coordinates x of the reference 
element {T,Tl}, one has to consider the chain rule and to 
implement 



3 (ur, o Y T , 



i(i) = J>W(i), 



'€{1,2,3}, (19) 



k=l 



(which follows from ( fT8] l) and 

d(T T ,), . . 2* . . d% 



d<Cj 



'■±(1) 



: E (a.) 

k=\ 



(f), {1,2,3}. (20) 



d20l ) leads to the Jacobian VT7; (i) of the geometric transfor- 
mation Tt, of ( fT~4-b . 

Considering material parameters such as Young's mod- 
ulus E and Poisson's ratio v, we can now compute with re- 
spect to the coordinates xeT the stress tensor, the von Mises 
stress (01, the elastic-plastic von Mises stress o v according to 
Neuber shakedown (O, the comparison strain £,, according to 
the Ramberg-Osgood equation (0|i and finally A^et according 
to the CMB approach (fTTb . 



1 

Ni7t 



dA = 



□ \^det(Y/;(s)) 



VgWjds (22) 



with charts y n (s) = T Ti (0,s h s 2 ),Ji2(s) = 
Y Ti ( 1 , * 1 , s 2 ) , Yi3 (s) r Ti (si,0,s 2 ), Y/4 (s) 

Y Ti (s u l,s 2 ) 7 yi 5 {s) = Tr.(5i,5 2 ,0),Yi6(s) = Y Ti (si,s 2 ,l) 
for s = (s\,S2) € [0,1] x [0,1] and with the Gram de- 
terminant g yi '(s) = det (VY/ / -(s) 7 "Vy^(s)), where Vy t j 
is the Jacobian of Jij. In case of tetrahedrons we inte- 
grate over the unit triangle A as map area with charts 
Yn(s) = r Ti (0,si,S2),ya.(s) = r Ti (s h 0,s 2 ),Yi3{s) = 
Yt;.(.si,.S2,0),Y/4(s) = Y Tj (s l ,s 2 ,l -Si -s 2 ). In case of 
volume integral the Gram determinant (|2H has to be replace 
by I det (DY7; ) I and the integration is conducted over the 
whole reference brick and tetrahedron, respectively. 

Now, we use quadrature formula (%i,(i>i)i=i /_ to nu- 
merically compute the cumulative distribution function (TT~2T >. 
Considering (|2TT i and (l22l leads to: 
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exp 



Net N F 
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n, 



■CO/ 



(23) 



where we defined the scale value 



'N d N F lq 

LEH^E (A / det(yi/fe))r 



-l/m 



CO/ 



(24) 




Fig. 3. FEA RESULTS OF ABAQUS 6.9-2 FOR THE VON MISES 
STRESS FIELD OF THE COMPRESSOR DISK. 

In case of bricks use Tableland in case of tetrahedrons Ta- 
ble 8.2 in Q0 . Equation ([23) shows that the Weibull haz- 
ard ansatz in the local and probabilistic LCF model leads 
to a Weibull distribution for the component Q. under surface 
driven and strain-controlled LCF failure mechanism. Here, 
the scale value r| is computed according to chosen finite el- 
ement types and quadrature formula. As mentioned in Sub- 
section fOl the parameters of the Ramberg-Osgood, the CMB 
equation and the shape parameter m have to be calibrated via 
LCF-test results. 

In this work the main input for FEA postprocessing are 
the coordinates and displacements of the nodes and the con- 
nectivity. Moreover, values for the parameters of our model 
must be given which depend on the considered material. 
Note that we do not employ additional information of FEA 
packages on which faces fn contribute to the surface 3£2 
with a non-vanishing area, so that we identify the surface 
only via the main FEA output. This identification of the sur- 
face can numerically fail if finite elements are significantly 
distorted. But in most applications there are at most a few of 
such distorted elements and their contribution to the overall 
probability of crack initiation can often be neglected. 



is subjected to an inhomogeneous stress field. The tempera- 
ture field is still homogeneous, but its value has increased. It 
is assumed that the fields are stationary which is only an ap- 
proximation regarding the real operating state of the disk. As 
a conservative estimation the temperature field of the shut- 
down state is set equal to that one of the final state, otherwise 
thermal mechanical fatigue (TMF) has to be considered. The 
transition from the shutdown state to the operating state and 
then back to the shutdown state is considered as one load cy- 
cle. It is assumed that the shutdown and operating state stay 
the same during the cycles which is again only an approxi- 
mation. 

An Abaqus FEA model has been created to predict the 
displacement and von Mises stress field of the compressor 
disk in the operating state. The model consists of approx- 
imately 10.000 Lagrange elements of the Serendipity class 
(C3D20) which is given in ( fT3l . The compressor disk is sub- 
jected to a homogeneous temperature field in the operating 
state which is considered in the FEA model by corresponding 
values for Young's modulus and Poisson's ratio. The shut- 
down state is already at our disposal as the displacement and 
stress fields are zero and the temperature field is set equal 
to the homogeneous one of the operating state. Thus, the 
main input for our postprocessor are the coordinates of the 
nodes, the displacements of the nodes and the connectivity 
of the FEA model prediction which determine the approxi- 
mation of the displacement field in the operating state and 
the geometry of the disk. 

Finally, values for the material parameters of the LCF 
model must be transferred. For this purpose we first cali- 
brate the LCF model with respect to the disk material ac- 
cording to standardized LCF tests. Considering the cumu- 
lative distribution function F^(n) of ( f23l and defining r| = 

(fdO.Nd™dA\ 1/,m , we obtain for the corresponding density 
function fy(n) the expression 



3 LCF-LIFE ESTIMATION FOR A COMPRESSOR 
DISK 

In this section we consider a linear elastic Abaqus model 
of a compressor disk, see Figure [5] and estimate its Weibull 
distribution with respect to LCF life by means of our post- 
processor. 

3.1 FEA Model and Parameter Calibration 

In the following we consider two states of the compres- 
sor disk: The shutdown state and the operating state. Before 
the gas turbine is activated the compressor disk is subjected 
to a homogeneous temperature field but to no stresses. This 
constitutes the shutdown state. In the operating state the disk 



d r \ m ( n\ m ( n\ m ~\ 

Mn) = Tn Un) = -{-) «*[-(-) j- (25) 



We subsume the parameters of the model in a vector 8, 
which includes the parameters of Ramberg-Osgood and of 
CMB and the Weibull shape parameter m. The experimen- 
tal data set for q strain-controlled LCF tests is given by 
{«,,£,-, 3i2,} (=1 . Here, n, is the number of cycles until 
crack initiation and e, is the strain on the gauge surface d£lj. 
We estimate (calibrate) by means of maximum likelihood. 



The log-likelihood function is defined as 

log(£({(n i ,e i ,3fl i )}te{i I ...,,})[e]) 

q (26) 
= £log(/ JV (3a i ,e i )(n i )[e]). 
(=1 

Let denote the likelihood estimator, then is given by 

Iog(£({(«,,Ei,3£l,)}fe{i , } )[0]) 

= max {log (L ({(n^dCl;)},^^) [9]) } • 

9 

Having optimized d27l i and found the estimator we can now 
apply our postprocessor to the FEA model of the compressor 
disk. 

3.2 Results of the Probabilistic Approach to LCF 

The most important result of the numerical integration in 
our model is the scale parameter r| of the Weibull distribution 
d23b which yields the probability for LCF crack initiation un- 
til cycle N. The corresponding Weibull shape parameter m is 
already estimated by the calibration of the previous section. 
The Weibull distribution is shown in Figure |4] for low num- 
bers of cycles compared to T|. Here, failure is defined by the 
initiation of the first LCF crack on the component and PoF 
denotes probability of failure and N* multiples of T|. 

For N* = 3.231 • 1(T 3 the PoF is 6.142 - 10~ 3 %, for ex- 
ample. If we consider that the compressor disk consists of 44 
disk segments, where one is shown in Figure [3] the probabil- 
ity for the initiation of the first LCF crack on the compres- 
sors disk is 0.270%. From a design perspective one decides 
which PoF is acceptable and then chooses the corresponding 
number of allowable shutdown and service cycles. Figure |5] 
shows the crack initiation density at N* on the disk's surface. 
This density corresponds to the local expectation value for 
the number of crack initiations. One can see that the crack 
initiation density is very much localized at the bearing flank. 
As we first compute the Weibull scale for every boundary 
face we can quantify this localization. In case of the previ- 
ous value for A^*, the 21 faces with the greatest crack initia- 
tion density have a combined PoF of already 5.531 • 10~ 3 % 
which is more than 90% of the PoF for all boundary faces. 

With respect to field data of already operating compres- 
sor disks no failure has occurred so far. Our postprocessor 
predicts for that number of disks and service cycles a low 
PoF. Note that the assumptions for the shutdown and final 
state as well as the FEA model only approximate the real 
operating conditions of these compressor disks. In particu- 
lar, consider that the field data consists of information from 



Probability of LCF Crack Initiation 




0.000 0.002 0.004 0.006 0.008 0.010 
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Fig. 4. WEIBULL DISTRIBUTION FOR LCF CRACK INITIATION 
ON THE COMPRESSOR DISK. 




Fig. 5. CRACK INITIATION DENSITY. 



different gas turbines whose service conditions can be very 
different. Moreover, uncertainties in the model parameters 

- recall the previous section and confer the end of this section 

- influence the real PoF of the compressor disk. Nevertheless 
the model is able to predict a low PoF. 

As a computational validation item of our tool we in- 
vestigated whether the numerical integration converges re- 
garding the order of the chosen quadrature. Since the FEA 
model consists of nonlinear Lagrange elements of Serendip- 
ity class the postprocessor employs quadratures of Table [T] 
confer Subsection l2.2l We computed the Weibull scale r| for 
every quadrature till the order of k = 11, i.e. till l q = 6 in- 



Convergence of Quadrature 




^ 1 1 1 1 

1 2 3 4 5 6 

Number of One-Dimensional Integration Points 

Fig. 6. CONVERGENCE OF WEIBULL SCALE REGARDING 
NUMBER OF INTEGRATION POINTS. 

tegration points in every dimension and thus 36 integration 
points on the unit rectangle. Figure [6] shows the results for 
each Weibull scale divided by the estimated limit and de- 
pending on the number i = 1 , . . . , 6 of integration points in 
single dimension. Using 16 integration points on the unit 
rectangle, approximately results in a converged value for r| 
which indicates a highly nonlinear behavior of the integrand 
N7l? of the probability function (l23l . This also justifies the 
choice of quadratures of higher order. 

The cumulative distribution function F^(n) in d23l de- 
pends on the parameters which are calibrated as described 
in Subsection 13.11 Because this calibration is a statistical 
estimation for depending on LCF test data, there are uncer- 
tainties for the values of 0. This effects the real PoF but is 
not considered by our method so far. Taking this additional 
uncertainty mathematically into account can be realized by 
Bootstrap methods or Maximum Likelihood asymptotic the- 
ory, for example. Confer iTPJI for more details on these meth- 
ods. 

Additional to the consideration of uncertainties in the 
model parameters 0, the extension of our model to inhomo- 
geneous temperature fields will be important for risk estima- 
tion for LCF crack initiation on components such as turbine 
blades. For this purpose a reliable temperature model for 
LCF is needed. Also note that the model could be extended 
to consider HCF, TMF and non-stationary FEA. From a de- 
sign perspective there is also the interesting possibility to op- 
timize the PoF ( 1231 with respect to the shape £2, i.e. find a 
design £2 under certain constraints such that the surface inte- 



gral and PoF d23l is minimized. This is also called optimal 
reliability, confer J5]. Because the integrand is sufficiently 
regular under additional smoothness assumptions there is 
even a link to gradient-based shape optimization, confer ll23l 
and l24l . This could accelerate computational optimization 
efforts significantly. 



Acknowledgements 

This work has been supported by the German federal 
ministry of economic affairs BMWi via an AG Turbo grant. 
We wish to thank the gasturbine technology department of 
the Siemens AG for stimulating discussions and many help- 
ful suggestions. 



References 

[1] M. Baker, H. Harders and J. Rosier, Mechanisches Ver- 
halten der Werkstoffe, third edition, Vieweg+Teubner, 
Wiesbaden, 2008. 

[2] D. Radaj and M. Vormwald, Ermiidungsfestigkeit, third 
edition, Springer, Berlin Heidelberg, 2007. 

[3] B. Fedelich, A stochastic theory for the problem of 
multiple surface crack coalescence, International Jour- 
nal of Fracture, 91, pp. 23-45, 1998. 

[4] D. Sornette, T. Magnin and Y. Brechet, The Physical 
Origin of the Coffin-Manson Law in Low-Cycle Fa- 
tigue, Europhys. Lett., 20 (5), pp. 433-438, 1992. 

[5] H. Gottschalk and S. Schmitz, Optimal Reliability in 
Design for Fatigue Life, Part I - Existence of Optimal 
Shaves, \sxXiv: 1210.49541 2012. 

[6] S. Schmitz, T. Seibel, T. Beck, R. Rollmann, R. Krause 
and H. Gottschalk, A Probabilistic Model For LCF, in 
preparation 

[7] P. Ciarlet, Mathematical Elasticity - Volume I: Three- 
Dimensional Elasticity, Studies in Mathematics and 
its Applications, Vol. 20, North-Holland, Amsterdam, 
1988 

[8] A. Ern and J.-L. Guermond, Theory and Practice of Fi- 
nite Elements, Springer, New York, 2004. 
[9] M. Knop, R. Jones, L. Molent, L. Wang, On Glinka 
and Neuber methods for calculating notch tip strains 
under cyclic load spectra, International Journal of Fa- 
tigue, Vol. 22, pp. 743-755, 2000. 

[10] M. Sherman, Spatial Statistics and Spatio-Temporal 
Data: Covariance Functions and Directional Properties, 
Wiley Series in Probability and Statistics, 2010. 

[11] A. Baddeley, P. Gregori, J. Mateu, R. Stoica and D. 
Stoyan, editors, Case Studies in Spatial Point Process 
Modeling, Lecture Notes in Statistics, 185, Springer, 
2006. 



[12] A. Klenke, Wahrscheinlichkeitstheorie, Springer, 
Berlin, 2008. 

[13] L. A. Escobar and W. Q. Meeker, Statistical Meth- 
ods for Reliability Data, Wiley-Interscience Publica- 
tion, New York, 1998. 

[14] P. Ciarlet, Basic Error Estimates for Elliptic Prob- 
lems, Vol. II: Finite Element Methods, ch. 2. Hand- 
book of Numerical Analysis, North-Holland, Amster- 
dam, 1991. 

[15] D. Braess, Finite Elemente - Theorie, schnelle Loser 
und Anwendungen in der Elastizitatstheorie, fourth edi- 
tion, Springer, Berlin, 2007. 

[16] R. B. Hetnarski and M. Reza Eslami, Thermal Stresses 

- Advanced Theory and Applications, Solid Mechanics 
and Its Applications, Vol. 158, Springer, Berlin, 2009. 

[17] W. Ramberg and W. R. Osgood, Description of Stress- 
Strain Curves by Three Parameters, Technical Notes 

- National Advisory Committee For Aeronautics, No. 
902, Washington DC, 1943 

[18] H. Neuber, Theory of Stress Concentration for Shear- 
Strained Prismatical Bodies with Arbitrary Nonlinear 
Stress-Strain Law, J. Appl. Mech. 26, 544, 1961. 

[19] M. Hoffmann and T. Seeger, A Generalized Method for 
Estimating Elastic-Plastic Notch Stresses and Strains, 
Part 1: Theory, Journal of Engineering Materials and 
Technology, 107, pp. 250-254, 1985. 

[20] G. Schott, Werkstoffermudung - Ermiidungsfestigkeit, 
Deutscher Verlag fur Grundstoffindustrie, forth edition, 
Stuttgart, 1997. 

[21] H.-O. Georgii, Stochastics - Introduction to Probability 
and Statistics, de Gruyter, Berlin, 2008. 

[22] O. Forster, Analysis 3, seventh edition, 
Vieweg+Teubner Verlag, Wiesbaden, 2012. 

[23] J. Haslinger and R. A. E. Makinen, Introduction to 
Shape Optimization - Theory, Approximation and Com- 
putation, SIAM - Advances in Design and Control, 
2003. 

[24] J. Sokolowski and J. -P. Zolesio, Introduction to Shape 
Optimization - Shape Sensivity Analysis, first edition. 
Springer, Berlin Heidelberg, 1992. 



